A1 = readmatrix('LD_Domain_1.txt')

x = A1(:,1)*2*pi/360
y = abs(A1(:,2))
z = A1(:,3)
y= y./z

x1 = A1(1:13,1)*2*pi/360
y1 = abs(A1(1:13,2))
z1 = A1(1:13,3)
y1= y1./z1

x2 = A1(13:36,1)*2*pi/360
y2 = abs(A1(13:36,2))
z2 = A1(13:36,3)
y2= y2./z2

x3 = A1(36:58,1)*2*pi/360
y3 = abs(A1(36:58,2))
z3 = A1(36:58,3)
y3= y3./z3

x4 = A1(58:80,1)*2*pi/360
y4 = abs(A1(58:80,2))
z4 = A1(58:80,3)
y4= y4./z4

x5 = A1(80:91,1)*2*pi/360
y5 = abs(A1(80:91,2))
z5 = A1(80:91,3)
y5= y5./z5

ax = polaraxes
polarplot(x1,y1,LineWidth=2,Color=[0 0 0.8])
hold on
polarplot(x2,y2,LineWidth=2,Color=[0.8 0 0])
polarplot(x3,y3,LineWidth=2,Color=[0 0 0.8])
polarplot(x4,y4,LineWidth=2,Color=[0.8 0 0])
polarplot(x5,y5,LineWidth=2,Color=[0 0 0.8])
polarscatter(x,y,5,'MarkerEdgeColor',[0 0 0])
%ax.ThetaZeroLocation = 'top'
%ax.ThetaDir = 'clockwise'
hold off

set(gcf,'Position',[0, 0, 150, 150])